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Abstract 

We present new results of existence of global solutions for a class of reaction cross-diffusion systems of 
two equations presenting a cross-diffusion term in the first equation, and possibly presenting a self-diffusion 
term in any (or both) of the two equations. This class of systems arises in Population Dynamics, and notably 
includes the triangular SKT system. In particular, we recover and extend existing results for the triangular 
SKT system. Our proof relies on entropy and duality methods. 


1 Introduction 

The purpose of this paper is to investigate existence and some properties of the solutions of the system 


OtU 

+ daU°‘ + dpv^) u]=u {Xy -XaU°' - Xb v’‘) 

in R+ X 17, 

(1) 

dtv - 

- Ax[{dy + dyV'^) u] = U (r„ - XyV^" - Xd u'^) 

in R+ X 17, 

(2) 


VxU ■ n = VxV ■ n = 0 

on R+ X 917, 

(3) 


u{0,-)=Uin, v{0,-) = Vin 

in 17, 

(4) 


where u = u(t,x) > 0, v = v(t,x) > 0 are the unknowns, the variables {t,x) browse R_|_ x O with O a bounded 
domain of R™ (m > 1), n = n{x) stands for the outward normal at point x of the boundary 917, Uin and Vin are 
nonnegative initial data, and the remaining terms are noimegative constant parameters satisfying 

D := {du,dy,da,df 3 ,d^,ru,rv,ra,rb,rc,rd,a,b,c,d,a,P,j} G x R+ x R^ x R+, 

(a > 0, d < 2 -h a, a < 1 + a) or (a = 0, d < 2, a < 1). 

The origin of this system is to be found in a well-studied system arising in Population Dynamics, known 
in the literature as the SKT system. The SKT system was introduced in [24] to model spatial segregation in 
two competing species of (let us say) animals (see also [21]). It has since then attracted the interest of many 
mathematicians, leading to a rich literature on the question of the existence of solutions (see for example [6] and 
references therein) and on the analysis of equilibria and stability (patterns are shown to appear; see for example 
[15]). Writing u and v the respective densities of the two different species, it takes the following form 

dtu — Aa;[(d„ + daU + dpv) u] = u (xy — XyU — r^v) in R_|_ x 17, 

dtv — Aj, [{dy + d^v + dsu) v] = v (xy — ryV — x^ u) in R_|_ x 17, (6) 

VxU ■ n = VxV ■ n = 0 on R+ x 917. 

When du = dy = da = dp = d-y = ds = 0, the system (6) reduces to the standard Lotka-Volterra competition 
ODS. The terms Xu, Xy are the intrinsic growth of the species, while Xb and Xd measure the demographic effect of 
the interspecific competition, and Xy, Xc indicate the demographic effect of the intraspecific competition. When 
non-zero, the terms Ax[{du + daU dpv)u\ and Ax[{dy + d-yV + dsu)v] model the spatial movements of the 
individuals in the domain 17. The positive constants du-, dy are standard diffusion rates, which indicate the 
frequency of the random walk of the individuals inside each species. The nonnegative constants dp and dg are 
usually referred to as "cross-diffusion" coefficients, and ecologically measure the repulsive effect, on the individuals 
of one species, of the presence of the individuals of the other species (as a result of the interspecific competitive 
pressure). The nonnegative constants da and dy, referred to as "self-diffusion" coefficients, measure the repulsive 
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effect on the individuals of the presence of the individuals of the same species (as a result of the intraspecific 
competitive pressure). 

The system (6) is often called "triangular" when ds = 0. From the point of view of modeling, this means that 
only the individuals of the first species tend to avoid the individuals of the other species. It therefore takes the 
form 


dtu - Ax[{du + daU + dpv) u\=u{ru-raU- n v) in IR+ x fl, 

dtv — Ax[{dv + d^v) v] = V (ry — TcV — Tdu) in K.+x fl, (7) 

VxU ■ n = VxV ■ n = 0 on IR+ x dil. 

In this case, the second equation is coupled to the first one only through zeroth-order terms (reaction), while in the 
full system (6) with ds > 0 both equations are coupled to the other one through both zeroth-order (reaction) and 
second-order terms (cross-diffusion). The full system (6) has a completely different structure from the triangular 
case (in particular it is possible to exhibit an entropy structure when ds > 0, see [6] and [11], but this entropy 
structure degenerates when ds = 0). 

In our study, we consider the class of systems (l)-(3). That is, we focus on a triangular type of cross-diffusion 
{ds = 0) as in the system (7), but in contrast to the system (7) where the diffusion rates {du + daU + dpv, dy -|- 
d-yV + dsu) and the growth rates (r„ — VaU — rsv, Vy — rcV — rd u) are required to be linear functions of u and v, 
in (l)-(3) we allow these functions to be more general power laws (with suitably chosen powers, (5)). Note that 
the class of systems we consider includes (7) (when a = l3 = j = a = b = c= d= 1). 

We now present our main mathematical result for this class of systems. 


1.1 Main Theorem 


We clarify the notion of weak solution we will use in the 

Definition 1.1. Let Vt be a smooth bounded domain o/K.™ (m € and let D € (]R1()^^ x IR+ x x K.+ . Let 
Uin ■= Uin{x) > 0 and Vm ■= Vin{x) > 0 6e two functions lying in L^{rt). 

A couple of functions {u,v) such that u := u{t,x) > 0 and v := v{t,x) > 0, and lying in x 

n) X LJ^^(R_|_ X LI) is a (global) weak solution of (l)-(4) if 

Va: [{du + daU°‘ + d/3W^) u] , [{dy + dyv'^) w] e L]^^{R+ X Lt) 


and, for all test functions tfi, 1(2 G C'g(IR+ x LI), we have the identities 

POO p p poc p poo p 

- / {dtipi)u- / fjl{0,-)Uin+ / / V xtpl ■ X [{du + dau"^ + dgv^) u] = / 'lfiu{ 

Jo Jq Jq Jo Jn Jo Jn 

pOO p p pOO p pOO p 

- / {dt'f’2)v- / 'llj2{0,-)Vin + / / Vx1p2 ■ yx[{dv + djV'^)v] = / / 1(2 V {Vy - Ty - Vd u‘^). 

Jo Jn Jn Jo Jn Jo Jo, 

Note that the assumptions on Uin, Vm, u, v, ipi, 1(2 ensure that all integrals in the two identities above are finite. 

Our main result is contained in the 


{Xu-Ta U°- -n v'’), 
( 8 ) 


Theorem 1.2. Let Ll be a smooth bounded domain 0 / R™ (m G N*^. Let the coefficients of system ( 1 ) - ( 2 ) 
satisfy condition (5). Consider initial data Um > 0, Vi„ > 0 such that Uin € 1^(0), Vm € L°°{Ll). 

Then, 

i) there exists u = u{t,x) > 0, v = v{t,x) > 0 such that {u,v) € T^^“(]R+ x Ll) x L^c{^+ X O) and {u, v) is a 
(global) weak solution of system ( 1 ) - (4) in the sense of Definition 1 . 1 . 
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Furthermore, this solution (u, v) satisfies for all T > 0 


sup V < max < sup u, 
[o,T]xn [ n 

rT 


lie 


f [ <C{n,T,u,„,v,n,D), 

Jo Jn 


/o Jn 


< C(p, D) (for alio < p < oo), 


[ u{T) + [ [ V,[(l + u)“/2] " <C(r!,T,Ui„,u,„,D) {tfa>0), 

JQ, Jo Jq, 

[ uiT)+ [ [ \\74\ogil + u)]f <C{n,T,u,r.,Vin,D) {ifa = 0), 

t/ Q J 0 J ^ 


(9) 

( 10 ) 

( 11 ) 

( 12 ) 

(13) 


where the constant C(0, T, D) only depends on the domain 12 (and the dimension m), the time T, 

the norms of the initial data ||Mm||L2(a) ||i^m||L°°(a) o,nd the choice of parameters D, and the constant 

C{p,fl,T,Uin,Vim D) only depends on the same quantities and the parameter p. 

Ifuin furthermore satisfies Uin{x) > 0 a.e. on 12 and logUm G L^(12), resp., if Vm furthermore satisfies Vin{x) > 0 
a.e. on 12 and \ogVin G L^(12), then for all T > 0 



[ 1 logMl(r) + du 

[ f |V,,[l0gM]|^ < 

/ |logUi„| 


(14) 


Jn .j 

r 

'0 Jn 
r^' r 

In 

r 



resp., 

/ |logu|(T) + dy 
Jn 

/ / |V,r[logu]| < 

Jo Jn 

/ |logn„| 
Jn 

+ C{9.,T,U^n,Vin, D). 

(15) 


ii) If a = 0, we have the additional estimate, for some v = i/(12,Ui„, D) > 0 depending only on the domain 12 (and 
m), the norm ||uin||Lcxj(Q^ and the parameters D, for all T > 0, and for some Ci (12, T, Um, D) > 0 depending 
only on the domain 12 (and m), the time T, the norms (||wi„||L 2 (Q), ||uin||L°°(n)) and the parameters D, 



^ C\ (12, T, Uin, Vin-, C^)- 


(16) 


iii) 7/7 = 0, assuming furthermore that Vm G W^’'^(12) for some q > 1 satisfying 


1 < q < {2 + a)/d if a > 0, 1 < q < {2 + i')/d if a = 0, 


(and assuming furthermore the compatibility condition \IxVin ■ n = 0 on 912 if q > 3), we have the additional 
estimate, for all T > 0, 

[ [ \dtv\^ + [ / iV't-r + / / |V,,u|'«<C'2(g,12,r,u„,u„,D), (17) 

Jo Jq Jo Jn Jo Jn 

where the constant C^iq, 12, T, Um, Vm, D) only depends on the parameter q, the domain 12 (and the dimension m), 
the time T, the norms of the initial data (||uiri||L 2 (Q), ||ui„||Lcxjp^ 2 , 9 (Q^) and the parameters D. 

We list in the remark below some possible extensions of the results mentioned in Theorem 1.2. 

Remark 1.3. Thanks to the bound (9) given by a maximum principle for equation (2), the power laws "v i—>■ " 

and "v I—>■ u**" in (1) can be replaced by any continuous function of v on K+ which is smooth (77^(1^^)^ and 
positive-valued on . The power law "v i—>■ v'^ " in (2) can be replaced by any continuous function on which 
is non-decreasing, smooth, and positive-valued ouM^. Following the same idea, but furthermore ensuring that the 
maximum principle remains valid for (2), the power law "v i— v‘^ " can be replaced by any continuous function on 
]R_i_ which is smooth, positive-valued on and which furthermore tends to +oo in +oo. With these replacements, 
all results of Theorem 1.2 hold, with the bound (9) adapted when necessary, and with the condition ^ = Q in iii) 
replaced by the condition that the function replacing "v i—u’*' " is constant. 

The power laws "u i—> " and "u i—>■ " can be replaced by any continuous function of u on ]R_|_, smooth and 
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positive-valued on and dominated by (or having the same behaviour as) "u i—>■ ", resp. "u i—> u‘^" in +oo. 
The power law "u " can be replaced by any non-decreasing continuous function of u on M+ which is smooth 

and positive-valued on and have the same behaviour as "u i—>■ " in +(X). With these replacements, all results 
of Theorem 1.2 hold. 

Assume a = 0. In this case, an estimate of type (10) ensures enough integrability (to get compactness) for the 
terms raU^'^°' in (1) and Vdvu'^ in (2) only when a < \ and d < 2. The slightly better estimate (16) is therefore 
crucial to estimate these terms in the two cases (a = 1, d < 2) and (a < 1, d = 2). It even treats these terms 
when a < l-\- u and d < 2-\-v. As a consequence, we expect the results of Theorem 1.2 to hold when the condition 
(a = Q, a < 1, d < 2) is replaced by the wider condition (a = 0, a < l-\- 1 ', d < 2 -\- v). Note that v can indeed be 
chosen independent of a and d (see Section 5.1). 

Assume a > 0 and 1 -\- a, d < 2 -\- a. As seen in Sections 2-f, in this case the proof of existence entirely relies 
on estimates of the type (9)-(13) (in particular we do not use estimate (16)J. As shown in the sequel, the proof 
of estimates of the type (9)-(13) only requires Uin to be in n H))j^(f2), where the space H~^(r2) is defined in 
Section 1.3. Therefore in the case a > 0 and \-\-a, d < 2-\-a, the results i) and Hi) for any 1 < q < {2-\-a)/d hold 
with the assumption Uin G L^(ll) replaced by the weaker assumption Uin G n H~^(f2) (and with the constants 
C and C 2 depending on ||um||LinH-i(n) instead 0/ 

Assume a = 0. In ii), estimate (16) is a consequence of a Lemma relying on duality techniques (namely. 
Lemma 5.1), which is adapted from a similar result from [5]. This Lemma can be somewhat improved. Indeed, 
following Remark 2.3 in [5], the constraint Um G L^(r2) can be replaced by the weaker constraint Uin G 
where vi is a small positive constant. Therefore, all results of Theorem 1.2 hold when Uin only belongs to the 
space (where vi only depends on Ll, m, ||i’in||L°°(f2) nnd D), with the constants Ci and C 2 depending on 

(f2) instead 0/ 

Assume 7 = 0. Estimate (17) being a direct consequence of the properties of the heat kernel (see Section 
5.2), we can replace the set W^’'^(0) in Hi) by the optimal set to apply the properties of the heat kernel, that is, 
the fractional Sobolev space Furthermore, the compatibility condition required for the case q = 3 is 

actually slightly weaker than 'N/xVm • n = 0 on dQ " (see for example [16]). 

1.2 In the literature 

In the last decades, mathematicians dedicated a considerable effort to the question of the existence of solutions 
for systems of the form (l)-(3), and particularly for the original system (7). 

The local (in time) existence of classical solutions was established by Amann in 1990 in the two papers [1], 
[2]. His theorem also provides a criterium to show that these solutions are global: it suffices to prove that the 
solutions do not blow-up in finite time in suitable Sobolev spaces. 

The global existence for the original system (7) has been investigated under various restrictive assumptions. 
Most results rely on Amann’s theory, therefore the problem is to prove bounds in appropriate Sobolev spaces. 
One of the main difficulties lies in the use of Sobolev embedding theorems in the parabolic estimates, which 
provide satisfactory results only in low dimension. Therefore, many existing results require strong restrictions on 
the dimension and/or on the parameters of the system (typically, one assumes that the cross-diffusion is weak, 
in the sense that the cross-diffusion term dp is small compared to some other parameters), see [7], [8], [17], [18], 
[19], [25], [26], [27], [28]. See [12] for a more detailed bibliography. 

So far, three groups managed to remove this type of assumptions on the dimension and/or parameters, in 
particular cases: Choi, Lui and Yamada in [8] (in the presence of self-diffusion in the first equation and in the 
absence of self-diffusion in the second equation, i.e. da > 0 and dj = 0) and, very recently, Desvillettes and 
Trescases in [12] (in the absence of self-diffusion in both equations, i.e. da = d-^ = 0), and Hoang, Nguyen and 
Phan in [14] (in the presence of self-diffusion in the first equation, i.e. da > 0). For the original system (7), our 
result is the first one treating the case {da = 0,dj > 0). Furthermore, it provides an unifying proof for the cases 
with and without self-diffusion (in one or both equations). 

We now mention some results of existence of global solutions for systems of form (l)-(3). Wang obtained it 
in [27] in the presence of self diffusion in the first equation {da > 0) and in the absence of self diffusion in the 
second equation {dj = 0), under a condition (depending on the dimension) of smallness of the parameter d w.r.t. 
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the parameter a. The case without self-diffusion {da = = 0) was solved by Pozio and Tesei in [23] under some 

strong assumption on the reaction coefficients, and by Yamada in [29] under the assumption a > d. We also 
mention the work of Murakawa [20] in which the reaction terms considered are Lipschitz continuous functions of 
u, V (and no self-diffusion appears). The results in the case without self-diffusion were extended by Desvillettes 
and Trescases in [12[. More precisely, in [12], the authors obtain global weak solutions for systems of the form 
(l)-(3) in the absence of self-diffusion {da = d-y = 0) with the following constraint on the parameters: (/3 > 1, 
a < 1, d < 2) or (/3 > 1, a < d). (Note that this indeed includes the original system (7) when da = dy = 0.) The 
main ingredients of the proof are entropy and duality methods. 

In the continuation of [12], the present paper deals with weak forms of solutions and exploits entropy and 
duality methods. These methods give rise to estimates for the solution which are quite explicit. Furthermore, 
considering weak forms of solutions allows us to consider initial data of low regularity (in comparison with most 
of previous works which rely on Amann’s theory). In comparison with [12], our Theorem includes the cases with 
self-diffusion terms {da > 0 and/or (i.y > 0). Another improvement is the possibility, for the first time, to consider 
cross-diffusion terms which are quite singular functions of v in zero, since we remove the assumption (3 > 1 
and replace it by the assumption f3 > 0. 

Finally, we refer to [11] for systems of the form (l)-(3) presenting in addition a cross-diffusion term in the 
second equation (non-triangular case). 

1.3 Notations 

When T > 0 is fixed, we write 17 = L^(]0,r[xn) for any p € [l,c)o]. For p S [l,oo[, we denote by = 
L^~''(]0,T[xn) the union of all spaces L*(]0,r[xn) for q > p. We write H/j(n) the space of functions of H^(r2) 
with zero-mean value on 11, and we write H~^(n) its dual space. 

We recall that D is the set of parameters, defined in (5). In the sequel, C{...) always denotes a positive 
constant, which only depends on its explicitly written arguments and may change from line to line. For example, 
C{n, D) is a positive constant that only depends on 12 and D. Furthermore, when it depends on D, the constant 
C{..., D) can be chosen continuous in D on the set {D € (K.;/)^® x R+ x K.^ x R+ : a<l +a, d<2 + a}. 

1.4 Plan of the paper 

We will first prove our result of existence under the following condition, which is slightly more restrictive than 

(5), 

a>0, d<2 + a, a<l + a. (18) 

More precisely, we will assume condition (18) in Section 4. 

In Section 2, we perform formal computations to establish a priori estimates on the solutions of system (1)- 
(4). In Section 3, we define a semi-discrete (in time) scheme designed to be a smooth approximation of system 
(l)-(4), and we prove (rigorously) estimates on the solution of this scheme that are independent of the time step. 
We use these estimates in Section 4 to pass to the limit in the approximating scheme under condition (18). This 
proves the existence of solution and the estimates required in i) when condition (18) is satisfied. Finally, we come 
to the end of the proof of Theorem 1.2 in Section 5. 

2 A priori estimates 

This section only contains formal computations. They will not be used in the proof of Theorem 1.2, but we 
believe that these computations are useful to clarify the main tools at the origin of our result of existence. Similar 
computations will be performed rigorously in Section 3 at the level of an approximating scheme. 

Therefore, suppose that {u,v) is a classical solution of system (l)-(4) for some smooth positive-valued initial 
data Uin = Uin{x) > 0 and Vin = Vin{x) > 0, and such that u and v are positive-valued and sufficiently smooth 
to perform all following computations (for example, in (7/(18.+ x fl)). Fix T > 0. Assume condition (5). In the 
following Sections 2.1-2.3, we compute estimates for the solution {u,v) on [0,T] x 17. 
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2.1 Basic estimates 

2.1.1 Maximum principle for v 

A direct application of the maximum principle to equation (2) gives 


sup[o,T] v{t) < max{supn Ui„, j }. 


(19) 


2.1.2 Mass conservation for u 

We integrate the equation for u on [0,T] x ft: 


rT 

/O J Vl J ^ JVt 

Since (for all u > 0) r^u < C{a,ru,ra) + we have 


[ u{T)= f Uin+ [ [ u{ru - raU°-- nv^). 

Jq Jq Jo Jq 


f^u{T) + ^fo < f^u^n + lfilTCXa,ru,ra) 


( 20 ) 


2.2 Duality estimate 

We now present an a priori estimate obtained thanks to a recent lemma relying on duality methods. This type 
of lemma was introduced in [22] and has since then showed to be very useful in the context of cross-diffusion (see 
[3], [4], [9], [10]). We cite here a version coming from [11]. 

Lemma 2.1. Let r„ > 0. Let M : [0, T] x O — >■ 1R+ be a positive continuous function lower bounded by a positive 
constant. Smooth nonnegative solutions of the differential inequality 

dtu — A{Mu) < r^u on fl, 
dn{Mu) = 0, on dLl, 


satisfy the bound 


where {um) denote the 

We check that M : 
equation (1) to get 


J J Mu^ < exp(2r„T) x ll^^mllH-i(n) A {um)^ J > 


mean value of x ^ u{0,x) on LI and Cq is the Poincare-Wirtinger constant. 

= du daU°‘ -|- dpv^ is lower bounded by du > 0. We therefore can apply the lemma to 



[du daU°^ dpyP] 




< C{n,T,Ur 


In particular, since all terms in the LHS are nonnegative, 



< C{n,T,u,n,ru) X ^1-t Jq J^[dudaU^^dfjvdfj 
< C{n,T,Uin,V,n,D) X (^1-t /p^/j.^]! -kdaU“]^ , 


where we used the L°° bound (19) in the last line. As a consequence, using furthermore the inequality (for £ > 0 
small enough, and for all z > 0) z°‘ < -\- £z'^~^°‘, we have 


lo In- C{Ll,T,Uin,Vin, D). 


( 21 ) 
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2.3 Entropy estimates 

We now present two new estimates which are crucial to obtain weak compactness on the solutions (as they yield 
a bound for the gradients of the solutions). These estimates are obtained thanks to the introduction of two 
functionals whose evolution along the flow of the solutions can be controlled. When decreasing, such functionals 
are called Lyapunov functionals and sometimes "entropies". By a small abuse, we will refer to the resulting 
estimates as "entropy estimates". 

2.3.1 Entropy estimate for v 

For any p ^ define 

E^(t) := / -{t). 

Jn P 

We compute the derivative 


dtEy{t) = [ dtvyP ^{t) 
[ vP{ry - rcv'^ - rdu‘^) + [ A^[{dy + d^v'*) v], 

Ja Jn 


where the last term writes 


[ yP ^A^[{dy+d^v'^)v] = -[Vx[v^ ^] ■ V^Udy + d^v'*) v] 
Jn Jn 

= - 4 ^/ \V,vP/^\^[idy+d^{j + l)v'>)]. 
P JQ, 


We integrate on f G [0,T]: 


/n P 


Now taking 0 < p < 1, we get 


-(T) + 4 


"-^1 ! 

P Jo Jq 


[idy+d-^{-f+l)v'^)] 


[ -( 0 ) + [ [ vPiry-ryV^-rdU^^). 

Jn P Jo Jn 


r [ V,K/ 2 ] " < C{p,D) I [ vP{T) + r [ + ) 

Jo Jn \Jn Jo Jn J 

and using the L°° bound (19) and the dual estimate (21) with d < 2 + a (given by condition (5)), 

[ f < C{p,n,T,Uin,Vin,D) (for all 0 < p < 1). 

Jo Jn 


( 22 ) 


Now, let g > 1. Let us pick some p g] 0, 1[. We can write = {q/p)v^'^ x G L°° x L'^ by (19) 

and (22). Therefore 


/o^/n < C{p,il,T,u^n,Vin,D) (for all 0 < p < cx)). 


(23) 


2.3.2 Entropy estimate for u 

Define 

We compute the derivative 


Eu{t) ■■= [ \og{l+u){t). 
Jn 
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dtE^{t) = [ ^{t) 

Jn 1 + M 

= / ^ (ru- Ta - n v'’) + / ——+ +(i/ 3 ti^)'u], 

in ^ in 

where the last term writes 

[ 7 -J— A^[((i„ + dc,u°‘ +di3V^)u] = - [ Va:[--^] • Vx[idu + daU°' + di3V>^)u] 

= [ \S/a:[\ogil + u)]\^ {du + da{l + a)u°‘+ df3V^)+ [ df 3 ---^V^[\og{l + u)]-y^v^. 
in Jq ^ + u 

We integrate on t G [0, T]: 

f log(l + M)(0) + [ [ \V4^og{l + u)]f {du + da{l + a)u°‘+di 3 V^) 

j d •i 0 'i o 

= /log(l + u)(r) - [ [ —^{ru-raU°^-nv^’) - [ [ dp—^V^[\og{l + u)]-V^v^. 

in io in 1 + u io in 1 + M 


Using the elementary inequality 2xy < , 
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/ / 1 , V4\og{l + u)] ■ V 

Jo Jn ^ + u 

= 2 

/ / V4l0g(l+u)]u^/2.v^^/3/2 

Jo Jq ^ + u 


1 

< - 

- 2 


/ / |V^[log(l + M)]|%^ + 2 / f 

/o in io in 


1 + M 




Reinserting in (24) (and using t6/(l + u) < 1), we get 


/ log(l + u)(0) + [ [ \V4^og{l + u)]f {du + dail + a)u°‘+ ^v^) 

J Q, J 0 Jq ^ 

rT 


< / log(l + M)(T) + 
in 


/ {ru + Ta + rt, u'') + 2dp / 
/o in io in 




(24) 


In the RHS, the first integral and the second term of the second integral are controlled thanks to the 
estimate (20) and the last term of the second integral is controlled thanks to the L°° bound (19) : 

/ log(l + m)(0 ) + [ [ \Vx[\og{l + u)]f {du + dail + a)u°‘+ 

J Q J 0 Jq ^ 

^ C Uin-, (1-, Tu-, Ta-, ‘2d^ 

Finally, the last term is controlled thanks to (23), so that 


r 


Jo Jn 



Io /a|V.[log(l+u)]|"(l + u“) <C(f2 ,T,Ui„,D). 


(25) 


3 Approximating scheme 

In this section we define a semi-discrete (in time) scheme intended to approximate system (l)-(4), and establish 
rigorously uniform estimates (w.r.t the time step) for the solution of this scheme. Thanks to these uniform 
estimates, we will be able to pass to the limit in the approximating scheme in the following Section (under 
condition (18)). 
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3.1 Definition of the scheme 

The scheme takes the following form: (uo,vo) are given and for 1 < k < N {N = T/t, T > 0) 

{{du da’ll dpVi^ )rr/i:] — Uk i^u '^k ^^^k) ^ ^5 


Vk — Vk-1 


- ls.^[{dv + d^vl)vk\ = Vk{rv-Tcvl-rdu'l), onVt, 
dnUk = dnVk = 0, on an. 


( 26 ) 

(27) 

(28) 


This scheme was introduced in [11] in a more general setting. More precisely, in [11] systems of the following 
form are studied : 


Uk — Uk-i 


- A 4 A{Uk)] = R{Uk), onn, 
dnA{Uk) = 0, on an, 


(29) 


where 


A : 


U = 


I 

■+ 

/ui 


\ui^ 


A{U) = 


/ai{U)u4 


\ai{U)ui^ 


and R 


satisfy the following assumptions: 

HI For all i, the functions Oj and are continuous from . 


U = 


to ^ 


I 

■+ 

/ui ' 


\ui/ 


RiU) = 



(30) 


H2 For all i, ai is lower bounded by some positive constant a > 0, and n is upper bounded by a positive 
constant r > 0. That is ai(U) > a > 0, ri(U) < r for all U € 

H3 A is a homeomorphism from 

Following [11], we introduce the 


to itself. 


Definition 3.1 (Strong solution). Assume HI, H2, H3. Let r > 0 a'nd let Uk-i be a nonnegative vector-valued 
function in L°°(n)^. A nonnegative vector-valued function Uk is a strong solution of (29) if Uk lies in L°°(n)^, 
A{Uk) lies in H^(n)^ and (29) is satisfied almost ever'ywhere on LI, resp. dfl. 

The general theorem from [11] writes 

Theorem 3.2. A ssume HI, H2, H3. Let LI be a bounded open set ofW^ 'with smooth boundary. Fix T > 0 
and an integer N large enough such that rr < 1/2, 'where r := T/N and r is the positive number defined in H2. 
Fix 7] > 0 and a vector-valued function L°°(n)^ 5 Uq > rj. Then there exists a sequence of positive vector-valued 
functions {Uk)i<k<N in L°°(n)^ which solve (29) (in the sense of Definition 3.1). Furthermore, it satisfies the 
following estimates : for all k > 1 and p G [1, oo[, 


UkG^^^y, 

Uk > Va.r.t on n, 

AiUk) G w 2 ’J’(n)^ 

where tja.r.t > 0 is a positive constant depending on the maps A and R and t, and 

/ (E^mI fe7l,(C/fc)) <Cin,Uo,A,r,NT), 


(31) 

(32) 

(33) 


( 34 ) 


fc=i 


^i=i 


where C{Ll,Uo, A,r, Nt) is a positive constant depending only on n, A, r, Nt and ]|C/o||LinH~i(n) • 
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( 35 ) 


Let us go back to system (26)-(28). It can be written in the form (29)-(30) with 1 = 2 and 

ai{u,v) = du + daU°‘ + dpv^, ri{u,v) = r„ - ra 

a 2 {u, v) = dv + djv'^, r 2 {u, v) = r„ — Vc — Vd u‘^. 

Applying Theorem 3.2 directly gives rise to the following existence theorem 

Corollary 3.3. Let Q be a bounded open set o/R™ with smooth boundary. Let the parameters of system (26)-(27) 
D G X R_|_ X R^ X R_|_. Fix T > 0 and an integer N large enough such that max(ru,r„)r < 1/2, where 

T := T/N. Fix rj > 0 and a couple of functions L“(n)^ 5 (uo,vo) > r/. Then there exists a sequence of couples of 
positive functions {uk,Vk)i<k<N m L°“(n)^ which solve (26)-(28) (in the sense of Definition 3.1). Furthermore, 
it satisfies the following estimates : for all k > 1 and p G [1, oo[, 

{uk,Vk) (36) 

Uk > ? 7 d.t, Vk > ? 7 d,t on Ft, (37) 

{[du + d^ul + dpvl,]uk, K + d^u/Jufc) G W2 >p(L!)2, (38) 

where po.r > 0 is a positive constant depending on D and r, and 

N 

V r / (1 + u^) < C{n, Uo, Vo, D, Nt), (39) 

fe=i 

where C{Fl,uo,vo, D,Nt) is a positive constant depending only on O, D, Nt and || (uq,' yo)||LinH-Tf 2 ) • 

Proof of Corollary 3.3. It suffices to check assumptions H1-H3 for A and i? defined by (35). Assumptions HI 
and H 2 are clearly satisfied with a = mm(du, dv) > 0 and r = max(rM,ri,) > 0. It remains to prove that the 
map A is a homeomorphism from R^ to itself. By a monotonicity argument, it is easy to see that the map A 
is a continuous bijection from R^ to itself. Since the map A is furthermore proper (thanks to the inequality 
|A('u,u)| > min((i„,(it,)|('u,u)| for all u,v > 0 ), it is a homeomorphism. □ 


3.2 Uniform estimates 

Let /i > 0 and let Uq = (uo,vo) be a couple of positive functions in L°°(n)^ such that uo,vo > /i > 0 a. e. on O. 
Let Lfk = {uk,Vk) be a solution of (26)-(28) (in the sense of Definition 3.1). We now establish uniform (w. r. t. 
N) estimates for Uk that will allow us to pass in the limit in the approximating scheme (in Section 4). Note that 
some of these estimates (namely, (40), (41), (54)) can be seen as the "discretized" (in time) version of the a priori 
estimates (19), (23), (25), while (39) can be seen as the "discretized" (in time) version of the a priori estimate 
(21). To establish rigorously these uniform estimates, we will use all the time the smoothness of Uk = {uk,Vk) 
(for any fixed N) specified by estimates (36)-(38). 


3.2.1 Maximum principle 
Lemma 3.4. We have 

max supufe < maxlsupuo, ( — I }. (40) 

k=i..N Q n \rc J 

Proof of Lemma 3.4- Recall equation (27). By the maximum principle, for all 1 < k < N, we have sup^ Vk < 
maxjsupQ Vk-i, ( ^ ) }) so that (40) follows by iteration on fc = 1..N. □ 


3.2.2 Entropy estimate for v 

Lemma 3.5. Assume condition (5). For all p G R+, we have 

N 

T 

fc=l 

|logUfc| + dv^T ( |Va;logUfe|^ < [ lloguol + C{FI,Nt,uo,vo,D), 
Jn jQ 


sup / 

0<k<N jQ 



Va-u 


,P/2 


< C{p,n,NT,uo,vo,D), 


(41) 

(42) 
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where the first constant only depends on p, O, Nt, ||Mo||LinH-i(n)) sup^ vq, D and the second constant only depends 
on n,NT, ||wo||LinH-i(n),supn vo, D. 

Proof of Lemma 3.5. Define for all p g] 0, 1[ and all 2 ; > 0 

zP 1 

(j)p{z) = z -1 + - > 0, 4>o{z) = z-log z > t). (43) 

p p 

We have (for all p s]0,1[ and all 2 ; > 0) 

4’'p{z) = 1 - zP-^, f'^iz) = 1 - 4>p{z) = {l-p)zP-‘^ > 0 , = ^ > 0- (44) 

For any p G [0,1[, by convexity of fp, for all p > 0, z > 0 we have (j)p{z){z — y) > (t>p{z) — (j)p{y). Multiplying 
equation (27) by 4>p{vk) and integrating over D, we therefore get 


/ [(ppivk) - fpivk-i)] - T / (l)p{vk)/iS.^[{dv + djvl)vk] < T / 4>'p{vk)vk{ry-rcvl-rduf). (45) 

Jo, JQ, JQ, 

Since (uk^Vk) satisfy the homogeneous Neumann boundary conditions (28), we can rewrite the second term as 


-T (j)p{vk)/iS.^[{dy+djvDvk]=T V^[(j)p{vk)] ■ Va;[id.u + d^vl)vk] 

JQ JQ, 

= T I fpMdylVxVkl"^ +T I (j)p{vk)dj{j + l)vl\VxVk\‘^. 
JQ JQ 

The last term being nonnegative, we have 

-T (j)p{vk)Ax[{dv + djvl)vk]>T (l)p{vk)dy\V xVk\‘^, 

Jn Jn 

SO that reinserting in (45) and summing for k = 1..N 


N N 

/ 4>p{vN) + '^r (j)p{vk)dy\'\/^Vkf < / 4>p{vo) + '^T (t>'p(vk)vk(ry - rcvl-ra 
JQ k—1 k—1 ^ 


-rd ui). 


(46) 

(47) 

(48) 

(49) 


Thanks to estimate (39) and the assumption d < 2 + a (consequence of condition (5)), we can control uj. in , 
so that, using furthermore the L°° bound (40) and the continuity of u >-)• vf^iv) on R+, we get 


/ (I>p{vn)+ '^'r (l)p{vk)dv\'^xVk\‘^ < / (l>p{vo) + C{p,fl,NT,Uo,Vo,D). 

Using the definition of fo and fp for 0 < p < 1, we have 

[ (j)p{Vk)dy\VxVk\‘^ = [ (I - p)vl~'^dy\VxVk\‘^ = [ {I -p)dy^\VxVl^‘^\‘^, 

Jn Jn Jn P 

/ (l>o{vk)dy\VxVk\‘^ = / v'f‘^dy\VxVk\'^ = / d„|Va;[logUfe]p, 
Jfi Jq, Jq, 


(50) 

(51) 

(52) 


so that (50) with 0 < p < 1 gives (41) for p < 1 and, using furthermore the elementary inequality (for all z > 0) 

I logz:| < z — logz, (50) with p = 0 gives (42). 

It remains to show (41) for p > 1. Let p > 1, and let us fix p g] 0, 1[. Combining (40) and (41) with p replaced 
by p, we have 




< C{p,n,NT,Uo,Vo, D). 


(53) 

□ 
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3.2.3 Entropy estimate for u 


Lemma 3.6. We have 


sup uk + 


0<k<N JQ 


^ x'^k 
1 + Mfc 


(1 + Mfc) < C(0, Nt, Mo, Mo, D), 


N 

sup / llogMfcl + du^^T / |Va;logMfe|^ < / | log Mo | + C {D,, Nt, Uq, Vq, D), 
o<k<N Jn Jq Jn 


where the constant C{fl, Nt,uo,vo, D) only depends on n,NT, ||Mo||LinH-i(n)) supfj mo, D. 

Proof of Lemma 3.6. Let (j){z) =2z — log(y^ + z) for all z > 0, with /r = 0 or /x = 1. It is useful to compute 

fiz) = 2z- log(^ + z), f'iz) = 2 - (j)'\z) = ^ > 0. (56) 

tk + z [fj, + zy 

By convexity of (j), for all y > 0, z > 0 we have (j)'{z){z — y) > (f>{z) — (f>{y). Multiplying equation (26) by (j>'{uk) 
and integrating over fl, we therefore get 

/ [^(Mfc) - ^(Mfc_i)] - T / (j)'{uk)/ls.x[{du + daU^ + di3v'^)uk] < T f' {uk)uk {tu - ra ul - rb v^) . (57) 

JQ, JQ, JQ, 

Thanks to the homogeneous Neumann boundary conditions (28), we can rewrite the second term as 

-T (j)'{uk)/lS.y{du+daU^ + di3v'^)uk]=T Vx(p'(uk) ' Vy{du + daU^ + djsv'^)uk] (58) 

Jfl JQ 

= T yx4>\uk) ■ [{du + da{l + a)ul + d/juf)Va;Mfe] +T Vxf)' {uk) ' [ukV xdjjvl] (59) 

JQ. Jo, 

= T (j)" {Uk){du + da{l + a)ul P dpvl)\\I xUk\^ + T / (j)"{uk)ukVxUk-yxdj3vl. (60) 

J Q JQ, 


Therefore 


/ [^(Mfc) - ^i'(Mfc-l)] + T / (l)"{uk){du +da{l +afu^ + dpv^)\VxUk\^ 

Jn Jn 

<-T J'iUk )Uk'\/xUk ■ xd/SV^ + T / <j)'{uk)Uk{ru - Taul - rbvl). 

Jn Jn 


The first term of the RHS can be rewritten 


-T f" (Uk)UkSJ xUk ■ ^ xdpv^ = <('"('Wfc)MfcVa;Mfc • d/3 2 (63) 

Jn Jn 

< ^ / <^"(Mfc)^Mfe|V 2 ,Mfcpd /3 uf + 2r / d /3 (64) 

z Jn Jn 

< ^ [ J'iuk)\yxUk\^d0V^ + 2T [ dp\Vxvy^f, (65) 

z Jn Jn 

where we used the elementary inequality 2ab < a^ + 6^ (for all a, 6 G R) and the bound 4)"{z) z^ = [z/(/i + z)]^ < 1 
(for all z > 0). Reinserting in (61), we get 

/" [(/)(Mfc) - (/)(Mfc_i)] + T /" (()"(Mfc)(d„ + da(l + a)Mfc + ^Mf)|Va;Mfcp (66) 

Jn Jn ^ 

<2t f d/s |Va;Mf^^P +T [ <l)'{uk)uk (r„ -VaU^- Vb v\). (67) 

JQ JQ 

Now using that (for all z > 0) (j)'{z)z{ru — ra z“) = (2 — l/(^ + z))z(r„ — To z“) < C{ru, Va, a), we have 

[ [(l>{uk) - (/)(Mfc_i)] + T [ (j)''{uk){du + da{l + a)ul + ^vl)\VxUk\^ (68) 

Jn Jn z 

<2t [ d /3 |Va,vf^^P + T|II|C'(r„,ra,a). (69) 

Jn 
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Summing for k = 1..N we get 


(I){un) + '^t [ (j)"{uk){du + da{l + a)Uk + )|Va;Mfcp 

' fe=i ^ 

r ^ r 

< / (j){uo) + 2 ^t / + NT\n\C{ru,ra,a), 

Jn JQ 


(70) 

(71) 


and using furthermore estimate (41) with p = P and p"{z) > 0 (for all z > 0), 


N 

[ (I){un)+'^t f <j)"{uk){du + dail + a)uk)\VxUk\'^ < ( (j){uo) + C{^1, Nt,uo,vo,D). (72) 

J O k—1 ^ ^ 

We conclude (54) by taking /r = 1 and using the inequality 2z > (j){z) = 2z — log(l + z) > z (for all 0 > 0), and 
we conclude (55) by taking ^ = 0 and using the inequality <j){z) =2z — \ogz > \ log 2 :| (for all z > 0). □ 


4 Proof of Existence 

We are now ready to pass to the limit in the approximating scheme (26)-(28) in order to obtain a solution of 
system (l)-(4), at least when the condition (18) is fulfilled. 

Proof of i) under condition (18). FixT > 0. Define the sequence (uo,Vo) = {uq ,Vq) := N) + l/N,Vin + 

1/fV), so that (uo,Vo) approximates (uin,Vi„) in (L^(D) n H~^(D)) x L°°(D) (when JV — >■ 00), and for all N, 
(uojVq) G L“(D)^ and (uo,vo) > l/N > 0. For any N large enough {N > 2T max(r„,r„)), use Corollary 3.3 to 
define the family of couple of functions {uk,Vk)i<k<N solving (26)-(28) with initial value {uq ,Vq). 

Note that (39), (40), (41), (54), (55), (42) are therefore valid with uq (resp. vq) replaced by Um (resp. Vm) in 
the constant C in the RHS. If logMi„ is in L^(r2), we furthermore have that loguo approximates logUm in L^(D), 
so that (55) actually yields a uniform bound (w. r. t. N). For the same reason, if logUi„ is in L^(D), then (42) 
actually yields a uniform bound (w. r. t. Af). 

Definition 4.1. For h := {hk)o<k<N o- given family of functions defined on il, we denote by the step in time 
function defined on M x D by 


N 

h (t, x) . ^o(^)l]—T,0] (0 4" 'y ^ ^fc(^)i-](fc —l)T,fcT] (^) 4” ^Ar(x)l]7Vr,+C3o[(^)' 

fc=l 

Note that by definition, for all p G [1,(X)[, we have 

IIZi^l|Lr(]o.T[xa) = / \hk{x)\Pdx 

\k=i 

With this notation, for any family h := {hk)o<k<N of distributions on we have 

N-l 

dth^ = ho{x) (g) S-rit) + (g Skrit) in x f})/ (73) 

k^O 

In particular, we can rewrite equations (26)-(27) as 

N 

dtU^ = uo{x) g d-r{t) + X! (^Ax[{du + + dl3V^)'^k] + {Vu - raU% - rbvl)uk^ g 4r(t), (74) 

k^l 

N 

dtv^ = vo{x) g 6-r{t) + ^r(^A^[{dy + d^vl)vk] + (r^ - g 4r(t)- (75) 



13 



We want to pass to the limit when TV = T/r —>■ oo in the two equations above. Estimates (39), (40), (41), (54) 
and (55), (42) can be respectively rewritten as 


f [ {u^fil + {u^r)<C{n,Uin,V^n,D,T), 

Jo Jn 

/ \ Ijc 

sup sup < maxjsupuiri, ( — ) }, 

[0,T] f2 n \''"c J 

for all 0 < p < oo, [ [ < C{p,n,T,u^n,Vtn,D), 

Jo Jo, 


sup uY + / 

[o,T] Jn Jo Jn 


V xU 


N 


l-\- U 


N 


{l + {yYT)<C{^,T,U^n,Vrn,0), 


sup 

[ llogu^l 

+ du 

f f Va;l0gM^^< 

[ |log«o^| 

“h Cy (r^, r, , D) , 

[0,T] J 

'n 

" 

'0 Jn 

In 


sup 

[ |logv^ 

+ dv 

f f l^xlogy^l" < 

[ llog^^^l 

+ C{9.,T,Uin,Vin, D) 

[0,T] 

Jn 


Jo Jn 

Jn 



(76) 

(77) 

(78) 

(79) 

(80) 
(81) 


Note that all bounds give rise to estimates which are uniform with respect to N (with the additional assumption 
that loguin, resp. logUin, lies in L^(ll) for (80), resp. (81)). As a consequence of (76)-(79), we have uniformly 
w.r.t. N 

yN G l2+“, € L“, V^yY € L\ G \Y. (82) 

The first bound is a direct consequence of (76), while the third bound is obtained by writing V^y^ = (1 + u^) x 
YJxIjY/Y + 2jY)) G X lY thanks to estimates (76) and (79). The second and last bounds are a direct 

consequence of (77) and (78) with p = 2. Using furthermore condition (18), (82) leads to 

{ru-raiy^r-n{v^Y)u^ Yv - rYy^Y - rYy^Y)v^ (83) 

{du + dYu^Y + dYY^)Yu.^ G {dx + dYY^YY'' G L“, (84) 

where 1 + y= {2 + a)/ max(l + a,d) > 1 - 

Rewriting system (26)-(27) as the following equalities (which hold in .^'(]0, T[xU)) 

-{m^ - SrU^) = l^x[{du + da{u^T + C^/3 (w^)'^)M^] + Yu “ ra(M^)“ - rb{vYY)2jY^ (85) 

r 

- SrV^} = AYYu + dj{y^Y)y^] + Yv - rdy^Y - rd{y^Y)2L^^ ( 86 ) 

r 

where Sr '■ u{t, x) i—>■ u{t — t,x), we finally get 

-{m^ - G l1(]0,T[,W“2’\U)), - Srv^} G Li(]0,r[, W^’^U)). (87) 

T T 

The uniform (w.r.t. N) bounds (82), (87) are sufficient to apply a discrete version of Aubin-Lions lemma (see for 
example [13]) to get the strong convergences (up to a subsequence) when N —S’ oo 

y^ —>■ M > 0 in L^, y^ —> u > 0 in L^. ( 88 ) 


Let us first check that these strong convergences allow us to pass to the limit in the uniform estimates (76)-(81) 
to get estimates (9)-(15). We can directly pass to the limit in estimate (77), and we use Fatou’s lemma to pass 
to the limit in estimate (76) and in the first terms in estimates (79), (80) and (81). To pass to the limit in (78), 
we first notice that YYY^'^ is uniformly bounded in lY'^ thanks to estimate (77), so that it converges weakly in 
lY. We conclude by using the weak lower semi-continuity of the L^(]0, T[, W^’^(U)) norm. 

To compute the remaining limit in (79), it is convenient to notice that (1 -|- m“) ~ (1 -|- u)“ (in the sense 
that there exist Ca, Ca > 0 such that for m > 0, we have Cq,(1 -I- m)“ < (1 -I- u“) < C'q(1 -I- m)“), so that 
(79) actually yields a uniform bound in for (1 +y^Y~'^ = A^joY |Va;(l if a > 0 (resp. 
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for |V^u^|2(l + = |V:,log(l + if « = 0). Since (1 + (rgsp. log(l + u^)) is uniformly 

bounded in thanks to estimate (76), it converges weakly in L^. Using the weak lower semi-continuity of the 
L^(]0, T[, W^’^(U)) norm, we get the desired bound. 

To compute the remaining limit in (80), we first observe that, thanks to Poincare-Wirtinger’s inequality. 




so that (80) yields a uniform bound for logu^ in L^. Together with the strong convergence ( 88 ), this implies 
that logM'^ converges towards logu weakly in L^, which allows us to pass to the limit in the second term of (80) 
thanks to the weak lower semi-continuity of the norm. The same arguments allow us to pass to the limit in 
(81). 

It remains to check that {u,v) is a solution of (l)-(4) in the sense of Definition 1.1. Clearly {u,v) lies in 
^max(i+o,d) ^ j^oo ^j^^nks to estimates (9), (10) and condition (18). Furthermore, using estimates (9)-(13), it is 
classical to validate the following computations (thanks to regularized approximations of (1 -|-m)“/^ and v^) when 
a > 0 


eD 


eD 


V,[da 2 + - V,(l + u)“/2 u“/(l -f u)“/2-i e 

a , 


6 L° 




eL^ 


eL^ em 


V^[{du+di3V^)u] = — {du+di3V^){l+u)'^ Va;(l -|-m)“/^ +dj3'\/xv'^ U e L\ 


and, (thanks to regularized approximations of log(l -I- u) and v^) when a = 0, 


eL°° eL^ eL^ gL^ 

Vx[{du+da+ dliV^)u]= — {du+da+ dpV^){l+u) Vx'^Og{l+u) +dl3VxV^ u G L\ 

Therefore, we also have that Vx [{du + daU^ + dpv^) u], Vx [{d^ + d^v'^) u] lie in L^. Let tpi, ^p 2 be two test 
functions in ([0,T[xD). Let us first extend 'ipi, ip2 on 1R+ x D by zero; then we extend ipi, ^/;2 on R x D in such 
a way that ip 2 be in Cg(K x D). Testing equation (74) with ipi gives 


- [ [ - f i’i{-T,-)uo 

J — oo J Q, J Q, 

^ p ^ f 

= ^'r/ Axipiihr) \{du + dau'^ + di3V^)uk\+^T ipi{kT)uk{ru - Vaul-nvl), 


k=l 


k=l 


(89) 

(90) 


where all terms are well debned (all integrands are integrable on their respective domains of integration). We 
rewrite this equation as 


- [ f - [ i’i{-T,-)uo 

J — OO J Q. J Q 

pOO p pOO p 

/ / Ax'ipi [{du + daiu^)^^ + di3{v^)^)u^] + / -tpi (ru-raiu^T 

J—ooJQ, J—ooJCl 


where 


N 


:=^tpiikT,x)l](^k-i)r,kT]it) —^ x)l]o,T]it) in L°°(R x D). 




(91) 

(92) 


(93) 
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Note that we also have 


N 


^ Ax'ilji{kT, x)l](k-i)T,kT]i't) —)• Axtpi{t,x)l]o,T]{t) in L°°(]R x fi). 


(94) 


fc=i 


Therefore, thanks to the uniform integrability (given by (78)-(79)) and the strong convergences (88), we have the 
convergences when N —> oo 


' —OOn 




N 


/O Jil 




{ru-ra{u^T -n {v^f) 


/o Jn 


V’i'w(r„ -raU°' - nv^), 


' —OO n 


A^-ij)^ [{du + daiy/^T + dp{y^Y)u^ 


lO JQ. 


Axipi [{du + daU°' + dpv^) u 


where the last term can be rewritten (since Va, [{du + daU°‘ + dpv^) u] lies in L^) 

f f A^ipi [{du + dau’^ + dpv^) u\ = - [ f Vxi’i ■ Vx [{du + daU°‘ + dpv^) 

Jo Jn Jo Jn 


(95) 

(96) 

(97) 

(98) 


It remains to treat the term coming from the initial data in (91). By the dominated convergence theorem and 
the definition of ug, 


- / 'tpi{-T,-)uo -- / 'l/^l{0,-)Ur 

Jn Jn 


(99) 


Replacing these four convergences in (91), we get that (8) is satisfied. A straightforward density argument allows 
us to replace ■i/'i in (8) by any test function in ([0, T[xO). Performing very similar computations for equation 
(75) with the test function il’ 2 , we can finally show that {u, v) is a local (in time) weak solution solution of (l)-(4) 
on [o,r]. 

Performing iteratively this proof on [0, 2T], [0, ST], ..., we can define {u, v) on K+ x fl in such a way that {u, v) 
is a (global) weak solution of (l)-(4) and it satisfies estimates (9)-(13) for any T > 0. □ 


5 Special cases 

This section is devoted to the end of the proof of Theorem 1.2, that is, the proof of i) when condition (18) is not 
fulfilled and the proof of ii) and iii). 

Recall condition (5). "Condition (18) is not fulfilled" exactly means that 

a = 0 and {(a = 1, d < 2) or (a < 1, d = 2)}. (100) 

In the subsection below, we will prove together i) and ii) under the wider condition 

a = 0 and (a < 1, d < 2), (101) 

which is the condition required in ii) in our Theorem. Note that we therefore freely have a second proof of i) in 
the case a = 0 and a < 1, d < 2. 

The last subsection is devoted to the proof of iii). 

5.1 The case a = 0 

The following Lemma is crucial to establish estimate (16). It is adapted from a similar result in [12] (one main 
difference being that the version stated below tackles weaker forms of solutions), which is itself adapted from an 
original result in [5]. 
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Lemma 5.1. We consider T > 0, a bounded regular open set of R'" (m G q > 1, a function R in 

L*(]0, T[xr2), and a function M in L°°{]0,T[xD,) satisfying 

R{t, x) < K and 0 < mo < M{t, x) < mi for t g]0, T[ and x G fl, (102) 


for some constants K > 0, mg, mi > 0. Then, one can find v g]0, 1[ depending only on D, and the constants mg, 
mi, such that for any initial datum Uin in any nonnegative very weak solution u € T[xr2) of 

the system 

dtu — Ax{Mu) = R < K in]0, T[xr2, 

u{0,x) = Uin{x) in fl, (103) 

Vx{Mu){t, x) ■ n{x) = 0 on ]0, T[x9fi, 

(in the sense that for all test functions if € Cf ([0, T[xO) satisfying Vxf’ ■ n = 0 on [0,T] x did, the equality 


f ( udtip + ( Uinip{0, ■) + [ [ MuAxif + [ ( ifR = 0 

Jg Jn JQ Jo JQ Jg Jn 


is verified), satisfies 


l|w||L2+,.(]Q_7.[xn) < Ct (||Mm||L2(n) + K) , 
where the constant Ct > 0 depends only on O, T and the constants mg, mi- 


(104) 

(105) 


Proof of Lemma 5.1. The proof relies on the study of the dual problem. More precisely since here we consider 
very weak types of solutions for (103), the analysis can be done rigorously only through (the dual problem of) 
a "regularized" version of (103) (that is, a family of systems with smooth data and from which we can recover 
asymptotically the original system (103) in some sense). The first step of the proof is to define a "regularized" 
version of (103) and study its dual problem (111). The second step is to ensure that any solution u (of the 
original problem) considered is the very limit (in some sense specified later) of the solutions of the "regularized" 
problem (note that this is a result of uniqueness for the original problem). Finally, the third step is to establish 
an estimate of the type (105) for the solution of the "regularized" system, so that (105) follows after passage to 
the limit. 


First step: regularization and dual problem. Let (p®)o<£<i be a family of mollifiers on R’"+^ such that 
for all 0 < £ < 1 , 

p"> 0 , p" G C“(R’”+i), suppp"cB„+i( 0 ,£), [ [ p" = l. (106) 

JR JR'" 

We extend M by (mg +mi)/2 on a layer of width 1 around ]0,r[xn, then extend it by 0 on R™+^, that is, 
M{t,x) = {mg + mi)/2 if 0 < dist((t,x),]0,r[xfl) < 1, M{t,x) = 0 if 1 < dist((t,x),]0,T[xn), (107) 
and we define the family := [p® * ^]\[g tJxsT (convolution in R™+^ then restriction on [0,T] x O). Therefore 

for all £ G]0,1[, G C°°([0,r] xH), mg<M^<mi, (108) 

for all p G [1, oo[, —> M in LP{]0,T[xH) when £—0. (109) 

We also define a family (uf„)e of smooth functions on il, which are identically zero on a (£-dependent) neighbor¬ 
hood of dn, and which approximates uin in L^(n). 

We are now ready to consider the following problem (which can be seen as a "regularized" version of (103)) 

{ dtu^ — Ax{M^u^) = Rin ]0, r[xn, 

u^( 0 ,x) = <„(x) in ( 110 ) 

yxiM^u^)(t, x) ■ n(x) = 0 on ]0, T[xdil, 

Note that since is C°°, we can always write Ax[M'^u^] = AxM"^ + 2VxM^ ■ VxU"^ + AxU^ in the sense 
of distributions (for a distribution). By classical results of the Theory of Linear Parabolic Equations (see for 
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example Theorem 9.1, together with the final sentence in paragraph 9, in Chapter IV in [16]), there exists a 
unique satisfying dtu’^, G L^(]0,T[xll) which solves the problem (110) in the strong sense. 

We now introduce the dual problem 

{ dtv"^ + = / in [0, T] x U, 

u®(T, x) = 0 in n, (111) 

VxV^{t, x) ■ n{x) = 0 on [0, T] x dV,, 


where / is any smooth function on [0, T] x 11. Since and / are smooth, this problem has a unique classical 
solution u® G (^“([OjT] x 11) (see for example Theorem 5.3, Chapter IV in [16]). We claim that there exists 
1^1 G]0, 1[ depending only on 12, mo, mi such that for all r G [2 — i^i, 2 + i^i] and all e g]0, 1[, 

l|Aa;U^||L''(]o.T[xn) < C{n,mo,mi,p)\\f\\L’-(]o^TixQ) and ||u^(0, OIIlV^) < C'(ll>■^o,ml,p,T)||/||L'■(]o.T[xn)• 

( 112 ) 

Lemma 2.2 (together with Remark 2.3) in [5] states that for any r g]1,2[, if (with the notations of [5]) 


C 


rriQ+n 


mi — mo 


< 1 , 


(113) 


then (112) is true. The proof of Lemma 2.2 (and Remark 2.3) never uses the fact r < 2, so we can use that 
(113) =^> (112) for any r > 1. It therefore suffices to check (113) for |r — 2[ small. This is done in the case 
r < 2 in the proof of Lemma 3.2 in [5[. For r > 2, following the ideas of the proof of Lemma 3.2, an appropriate 
interpolation leads to the bound for all m > 0, 4 > r > 2 


Cm,r 


— ^ '^m,i 1 


, e 1-6 1 

where - + ^ 


that is, for all m > 0, 4 > r > 2, 


Therefore for all 2 < r < 4, 




(114) 

(115) 


^ mi - mo ^ 

O mp + mi --- <v 

2 I 


mo + mi 


1—4/r 


c 


,2-4/r 


mi — mo 


(116) 


The RHS is a numerical function depending only on f2, toq, TOi and r, and it tends to < 1 when r > 2 

tends to 2. Therefore there exists a small ixi > 0 depending only on 12, mo, mi such that the RHS is < 1 for all 
2 < r < 2 + ^ 1 . This implies (113) for all 2 < r < 2 + ixi. 


Second step: uniqueness. Let p G [2,2 + ^i[ and let m be a very weak solution of (103) such that u G L*’ 
(where l/p+ 1/p' = 1). We claim that u is the (unique) strong limit in of the family (M®)e defined in 

( 110 ); 

Since u® is a smooth function satisfying the homogeneous Neumann boundary condition, we can use it as a 
test function in (104) and against (110). Subtracting the two equalities thus obtained, we get 


f f {u^-u)dtv^+ [ {ul^-Uin)v%0,-)+ [ f {M^u^ - Mu)A,^v^ = 0. (117) 

Jo Jq Jq Jo Jn 

Using the smoothness of / and using u^,u G L^, we can perform rigorously the following computations: we 
use the definition of in (111) 

f f {u^ - u)dtv^ = - f f M^{u^ - u)AxV^ + f f{u^-u)f, 

Jo Jn Jo Jn Jo Jn 

and replace in (117) to get 

eT 


[ [{u^-u)f+ [ - M*„)u^(0, •) + [ [ {M^ - M)uAxV^ = 0. 

Jo Jn Jn Jo Jn 
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Using (112) with r = 2 + and Holder’s inequality, it yields 


r f iu^-u)f 

< 

/ Kn-Utn)v% 0 ,-) 

+ 

[ [ {M^-M)uA^v^ 

Jo Jn 


Jn 


Jo Jn 


< (7(12,mo,nil, T) + ||(Af -^)'w||L(2+i/i)'Qg ll/llL^+''i(]o,r[xO)' 

By duality, we obtain 

11"*^ ~ (]o,t[xO) — ^*(12, mo, nil, T) — UinWh^^Q) + ||(2U — M)u\\-^( 2 +vi)' (]o,T[xn)) ■ 

Using again Holder’s inequality, we end up with 

11^ ~ '^IIl(2+>'i)'(] o,T[xn) — C'(12, mo, mi, T) — Uin\\ij‘^(Q,) + ||2U — Af ||LS(-]o,T[xn)II'^^II lp'(] o,r[xO)) > 

where 1/(2 + vi) + 1/s = 1/p. Letting e tend to zero, we get that ||u^ — u||l( 2 +,.i)'(]q t’JxO) ^ 


Third step: uniform bound for the regularized problem. Let us choose / < 0 in (111). By the maximum 
principle, we have > 0. Since is a smooth function, we can use it as a test function against (110). We get 








As before, using the smoothness of u®, / and using S L^, we can rigorously compute 


(118) 







(119) 


and replace in (118) to get 


f [ u^f= [ <„u^(0,-)+ [ [ Rv^ < [ <„u®(0,-) + 2L [ [ v^, (120) 


where we have used the bound R < K and the nonnegativity of u®. The last term can be rewritten 


K r [ v^=K r f [ dtv^=K r f [ [-m^a,v ^+/]. 

Jo JQ Jo Jo Jn Jo Jo JQ 

Replacing in (120) and using (112) with r = 2 — and Holder’s inequality, it gives 

[ [ u^(-/) < C'(f^,wo,mi,T) (||M)’„||L2(fi) +iL) 

Jo Jn 


( 121 ) 


( 122 ) 


f 

Jo Jn 


Using the strong convergence of in and the smoothness of / to pass to the limit in the LHS, we get 

u{-f) < (7(12, mo,mi, T) (||Mm||L2(n) + K) ||/||l=— i(]o.T[xa)- (123) 

I 

Since u > 0 by assumption, this yields by duality 

(]o,T[xO) A (7(12,mo,mi,T) (||Mi„||L 2 (-f 2 ^ + at) . (124) 

This concludes the proof for some 0 < 12 < min{(2 — 1 ^ 1 )' — 2, 2 — (2 + 1 ^ 1 )'}. □ 

We are now ready to perform the 
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End of proof of i) and proof of ii). Let a = 0 and a, d satisfy condition (101). Let (oe)e C [a/2, a] and {de)e C 
[d/2,d] be two (strictly) increasing families of real numbers such that Ue —>■ a and de ^ d when e ^ 0. This 
implies that for all e > 0 , Oe < 1 = 1 + a and de < 2 = 2 + a, so that condition (18) is satisfied with (a,(i) 
replaced by (aejdg). Therefore, we can apply the results of Section 4 with this choice of parameters. For all 
e > 0, there exists a weak solution (ug > 0,V£ > 0) of (l)-(4) with {a,d) replaced by {agjdg). Furthermore, for 
any hxed T > 0, this solution satisfies estimates (9)-(13) (and (14), resp. (15) when logMi„, resp. logUi„, is in 
L^(0)) with {u,v) replaced by (ue,Ve). Since the constant C'(..., D) is chosen to be continuous in D on the set 
{D :a<l + Q;, (i<2 + a}, the estimates actually give uniform bounds w.r.t. e. As a consequence, we have the 
following uniform (w.r.t. e) bounds (for all p > 0 ) 

UeGL“, V,,log(l + ae) GL^, V^yPeL"^. (125) 

Let us check that we can apply Lemma 5.1 to Ug. We define Rg := (r„ —raw/' —rbVg)ug and Mg := du + da + dpv^. 
For all e > 0, i?£ G L'^' where Qg := 2/(1 + Og) > 1, and 


Rg < sup(r„ - raU°‘^)u 

u>0 


Qg / Tu 
Og + 1 V^o(ag + 1) 


1/Oe 



<K, 


where AT := r„ I ^) if > ra, AT := I ^ I if r„ < ra- Thanks to the bound (9) for Vg, we also have for 
all e > 0, Mg G A°° and 


0 < Too := du + da < Mg < d^ + da + dp max 


supr;f„, 

a 



=: nil < oo. 


We can therefore apply Lemma 5.1, which yields the bound for all £ > 0 


||ug||L 2 +-< (^(T, n,mo,mi) {\\uin\\L^(n) + K), (126) 

where v = ^(11, toq, rrii) > 0. Note that the constants K, niQ and toi being independent of e, this bound is also 
independent of e. 

From (125), (126), we have the uniform (w. r. t. £) bounds 


V^ug G G (127) 

{ru-rauf^-rbv'^)ug (r„ - CcU/- G (128) 

(du + da + dpv^)ug G l2+", {d. + djv])vg G L“, (129) 

where we used the computation V^Ug = (1 + Ug)Vx log(l + Ug) G x L^ for (127) and the assumptions on the 
parameters Og < a < 1 and dg < d < 2 for (128). Using the equations of {ug,Vg), estimates (128)-(129) yield a 
uniform (w. r. t. e) bound for dtUg, dtVg in L^“''^/^(]0, T[, W“^’^~'’'^^^(U)). Combined with the gradient estimates 
(127), this allows us to apply Aubin-Lions Lemma, so that, up to a subsequence, when e — 5 - 0, 

Ug —)• u > 0 in L^, Vg —^ u > 0 in L^. (130) 


We obtain estimates (9)-(15) for (u, v) with the same arguments as in the passage to the limit in estimates (76)- 
(81) in Section 4, and we obtain (16) thanks to (126) and Fatou’s lemma. As in Section 4, estimates (9)-(13) 
enable us to check that {u, v) G x L°° with Va;[(d„ + d^ + dpvd)u\ and Vx[{dv + d-yv'^)v] in L^. Using 

furthermore estimates (128) and (129), it is classical to check that {u,v) is a global weak solution of (l)-(4). □ 


5.2 The case 7 = 0 

Proof of in). When 7 = 0, the system satisfied (in the weak sense) by v can be rewritten as 

dtv — d'xAxV = f in K_|_ x U, (131) 

\7xV • n = 0 on R+ x dfl, (132) 

u(0, •) = Vin in O, (133) 

where d/ = d„ + d.y > 0, / = u (r„ — rgv'^ — r^u'^) G L/q^(K_|_ x U) (thanks to the estimates (9), (10) and (16)) 
and Vin G 1U^’*(U) (by assumption). Using the properties of the heat kernel, we get the L/,^(IR_|_ x U) bounds for 
dtV and stated in (17). The bound for VxV is obtained by interpolating (9) and the bound for V^u. □ 
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